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An accurate numerical consideration is carried out of the 
ground state for the simplified model which is traditionally 
used for the description of Verwey transition and related phe- 
nomena. In the framework of ID spinless fermion model, the 
effects of next-nearest neighbour (NNN) interaction on the 
metal-insulator transition are investigated for electron con- 
centrations 1/2 and 2/3. It is shown that for large enough 
NNN transfer integrals the electronic topological transition of 
metal-metal type is also possible. The corresponding phase 
diagrams are presented. 

PACS numbers: 71.30.-Fh, 71.27.-Fa, 71.20.Hk 



I. INTRODUCTION 

Studies of highly correlated electron systems are one 
of the "hot" areas of condensed matter physics. A 
tremendous number of papers have been concerned with 
the problems such as heavy fermions, Mott transitions, 
magnetism of highly correlated systems, unconventional 
mechanisms of superconductivity etc. In this context 
the classical problem of charge ordering which has been 
raised by Wigner 0] and Shubin and Vonsovsky |^ (see 
also 1^) received far less attention. As early as the 1930s, 
the famous Verwey transition in magnetite Fe304 at the 
temperature Ty = 120K has been discovered. According 
to the traditional views, this is the transition from the 
low-temperature insulating phase with Fe'^+ and Fe^+ 
ions forming superlattice (Wigner crystal) to the high- 
temperature conducting phase with the destroyed charge 
ordering Q. However, recent experimental data point to 
a more complicated physical picture since the long-range 
charge ordering below Ty appeared to be pronounced 
rather weakly [^Ug] and the infrared spectra show "tun- 
neling modes" which are presumably connected with the 
"smearing" of Fe'^"'" and Fe^+ ions over finite clusters M. 
The assertion that namely the charge ordering of Fe^ 
and Fe^+ ions is the cause of the metal-insulator transi- 
tion is also called in question in a number of experimental 
works d,^. Phenomena similar to the Verwey transi- 
tion are observed also in a number of other compounds, 
e.g., in R3X4 (R=Sm,Eu; X=S,Se), Yb4As3 (see |lO) and 
references therein), in the layered compounds RFe2 04 
,where R is a rare-earth ion Some of these com- 

pounds (e.g., EU3S4) exhibit long-range charge ordering 



at low temperatures, and other ones (e.g., Sm3Se4) show 
only short-range ordering. In the latter number 
of physical properties appeared to be anomalous, e.g., a 
giant low-temperature heat capacity which is not con- 
nected with current carriers is observed |ljjl^. Irkhin 
and Katsnelson jl^] suggested that an RVB (or pseu- 
dospin liquid) -type state with separated pseudospin and 
current degrees of freedom may be formed in such sys- 
tems. Hence, the problem of description of the Verwey 
transition and related phenomena seems to be rather in- 
teresting from both experimental and theoretical points 
of view and may turn out to be closely related to the 
more popular areas of the physics of highly correlated 
electron systems. 



II. THE FORMULATION OF THE MODEL 

To consider the effects of strong correlations in crystals 
including charge ordering the polar model ||^ is usually 
used with the Hamiltonian 



H = 



1 



ij JliTlj 



(1) 



where and Cio- are the creation and annihilation elec- 
tron operators, i is the site (Wannier state) index, a =t , | 
is the spin projection; = cf^Cia', Ui — + is the 
operator of electron number on site i; Uj is the transfer in- 
tegral; U, Vij are the matrix elements of Coulomb interac- 
tion on-site and between different sites correspondingly, 
and the primed summation symbol indicates that i ^ j. 
In comparison with the original polar model this expres- 
sion does not contain the multielectron contributions to 
the transfer processes and the direct exchange interaction 
(see ||] ). The effects of electron-phonon interaction may 
be partially taken into account by the renormalization of 
the model parameters (e.g. the formation of small-radius 
bipolarons by the condition V <Q etc.). 

In the case of Verwey transition in the magnetite the 
excess electrons distinguishing a Fe^+ site from a Fe'^"'" 
one are almost completely spin-polarized {Ty « Tc 
where Tc — 850i^ is the Curie temperature) and fill the 
band by less than half. Therefore the term with a ~[ 
may be omitted in the Hamiltonian (1) and one may pass 
to the so-called spinless fermion model 
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H : 



1 

- - 2^ V^jninj. 



(2) 



In the limit U ^ oo the Hamiltonian (1) coincides 
with (2) More exactly, this is valid only in ID 

case since in 2D and 3D systems the transfer Hamiltonian 
is renormalized by the spin-polaron effects |l6| . But if we 
deal with ferromagnets (this is the magnetite case) the 
latter are absent. 

The present work is devoted to the rigorous numeri- 
cal investigation of the ground state of the model (2) in 
ID case with the allowing for the interaction and trans- 
fer processes for the next-nearest neighbors. As it will 
be shown below they lead to a number of novel effects 
in particular to the opportunity of dimerization or elec- 
tronic phase transitions of " metal- metal" type. ID model 
may be considered not only as a simplified way to under- 
stand some features of real 3D case but also as a model 
to describe the processes in the conducti ng p olymers and 
other quasi-one-dimensional crystals |[l7|-p^. 

At tij — the model (2) is equivalent to the Ising 
model and describes in the ground state a "frozen" elec- 
tron distribution (Wigner crystal) depending on the spe- 
cific form of Vij. At Vij = we have the case of non- 
interacting electrons with the homogeneous distribution. 
Intuitive considerations of energy gain for competing pro- 
cesses show that at Uj ^ Vij the ground state should 
change. 

In the nearest-neighbor approximation ID model with 
the Hamiltonian (2) appears to be exactly soluble pO| , pl| . 
For the electron concentration per site p — 1/2 and V = 2 
(here and further we put ti,i±i — 1, ti^i±2 = t' , Vi^i±i — 
y, Vi,i±2 = V) the metal-insulator transition takes place 
with the appearance of the energy gap in the electron 
spectrum at V > 2. Let us investigate for the beginning 
the effect of the interaction V on this transition at t' — 0. 

For t = and p = 1/2 the ground state in the model 
(2) is determined by the ratio of V and V and, as it 
may be easily demonstrated, corresponds to the " usual" 
charge ordering [...10101010...] for V < and aimer 
ordering [...110011001100...] for V > \V (1 means the 
electron and means the hole on a site). So the model (2) 
is frustrated because of the competing nature of the in- 
teractions V and V (they favor different types of charge 
ordering) . 

Let us define the quantity Aoo as follows: it is equal 
to if the ground state degeneracy g is macroscopically 
large (i.e. increase faster than N when N — > cx)) at t = 0; 
otherwise it is equal to the difference between the energies 
of the lowest excited state and the ground one at i = 0. 
It may be shown by the direct enumeration of possible 
states with the same energy that at i = the degeneracy 
g is exponentially large iniV at p ~ 1/2, V^' = \V and 
therefore Aqo = in this case. 

One may obtain for p =1/2 



V - 2V' if \V > V, 



2V' ~V ii\V <V' <V, 



(3) 



I V 




if V > V. 

V iiV > V, 

V iiV < v. 



(4) 



As the electron transfer is taken into account the 
ground state for the case Aoo = is split into the energy 
band with the bandwidth of order of Nt. Since g » N 
this band should be continuous in the limit TV — > oo. 
Therefore one should expect that no energy gap will ap- 
pear between the ground state and the lowest excited one 
so the ground state will be always conducting at Aoo = 0. 
At Aoo > and t > it is naturally to expect that the 
conducting ground state will transform to the insulating 
one with the increasing Aoo- The following question is 
to be elucidated: at what value of A^(F, 1^') does the 
metal-insulator transition takes place and does the tran- 
sition point depend on V and V' separately? 

Hence, the following is the list of questions to be an- 
swered: (i) Does the metal-insulator transition exist at 
Aoo = 0? (ii) At p = 1/2, does the metal-insulator tran- 
sition occur at Aoo ~ 2 similarly to the exactly solu- 
ble case V' = 0? (iii) More generally, do the interaction 
constants enter the expression for the transition point 
mainly via the parameter Aoo? (iv) How does the metal- 
insulator transition proceed at p ^ 1/2 in particular at 
/9 = 2/3 (this filling corresponds to the situation in many 
real compounds, see e.g. [0)? 

These problems were considered already in a num- 
ber of papers (although, to our knowledge, for the case 
p = 1/2 only). In [g2[ the phase diagram of the ID 
spinless fermion model (rewritten in terms of spin op- 
erators po| ) with next-nearest-neighbor interaction has 
been obtained by the renormalization group approach. 
In [^,|^ the metal and insulator regions of this model 
has been investigated by an exact numerical calculation 
of " Drude" contribution to the frequency dependent con- 
ductivity (rather general approach to the numerical cal- 
culation of correlation functions for ID many-electron 
systems see e.g. also in [^^). However, the renor- 
malization group method can give only general shape of 
the phase diagram and not the exact numerical results for 
the positions of phase boundaries. On the other hand, 
the calculations in [Q have been carried out only for 
few sets of parameters of the model. The recent work 
contains the most detailed numerical investigation 
of the metal-insulator phase diagram in the model un- 
der consideration. Below, the problem of metal-insulator 
transition in the model (2) will be investigated by the 
calculations of the total energy and energy gap by the 
Lanczos method for the broader range of the parame- 
ters than in the previous papers and the answers to the 
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above questions will be sought for. Also the interesting 
phenomenon of the metal-metal transition which appears 
to be possible in the model with next-nearest-neighbor 
transfer processes will be studied. 

III. METAL-INSULATOR TRANSITION: THE 
RESULTS OF CALCULATIONS 

To investigate the ground state of the system we have 
used the Lanczos method of exact diagonalization for fi- 
nite clusters [^,^ (really we made the calculations for 
TV < 20) with the extrapolation to the limit N oo. 
The comparison of the results with the exact ones at 
t' = 0,V' = ||2^,0 have shown that we have at least 
four accurate digits in the ground state energy and two 
ones in the energy gap (which was calculated as the av- 
erage difference of the ground state energy at the adding 
and removing the electron) . Such accuracy is sufficient to 
answer the questions formulated above. The main results 
are the following: 

1. At Aoo = the metal-insulator transition is absent 
and the ground state is always metallic both for p = 1/2 
and for p = 2/3 (the calculations have been carried out 
up to 1/ - 102). 

2. The metal-insulator phase diagram for p = 1/2 
is shown in Fig.l. The metal-insulator transition takes 
place at A^o ~ 2. Qualitatively this result seems rather 
natural. Nevertheless, it is not trivial, to our opinion, 
that it turned out to be valid with the accuracy com- 
parable to the accuracy of the calculations. This result 
may be used for the testing of different approximations 
proposed for the investigation of 2D and 3D systems. 
Qualitatively our phase diagram is in agreement with 
that from . Moreover, the part of the phase diagram 
for which real calculations has been carried out in ||2j] 
is the same. But our results for the region of large V' 
are new. The main difference is that, according to our 
results, the boundary lines of metal-insulator transition 
are not crossed and the metallic phase is continued to the 
infinity along the line V' = V/2 in correspondence with 
the qualitative discussion in the Introduction. 

3. Fig. 2 shows the phase diagram for p = 2/3. The 
criterion of the metal-insulator transition is Aqo ~ 3. 
This result is new even for the case V = 0. The phase 
diagram for this case is qualitatively different from that 
for p = 1/2 (compare Figs.l and 2). In particular, the 
ground state is metallic for small enough V' and arbi- 
trary V. It is connected with the macroscopically large 
degeneracy of the ground state at V = 0. Also, it is not 
trivial that the position of the metal-insulator boundary 
with the accuracy of the calculations depend only on V 
at V > V and only on V at V > V . Since the ID 
spinless fermion model for p — 2/i was probably not 
considered earlier we present in the Table 1 the results 
for the ground state energy and the energy gap. 

4. We also have carried out the calculations for the 
case p = 1/2 with the adding of the third-neighbor in- 



teraction term y^^n^ni+a. It appears that the crite- 
rion of the metal-insulator transition Aoo ~ 2 is valid 
also in this case provided that the interaction potential 
is downwards-convex namely V + V" > 2V', Aoo being 
equal to 

{V-2V' + 2V" if \{V + V") > V\ 
2V' -V if i(l/ + V") <V' <V, (5) 

V if V > V, 

This statement is illustrated by the data of the Table 
2: in the case of the downwards- convex potentials (upper 
part of the Table) the criterion Aoo ~ 2 holds but no 
simple criterion can be established for the opposite case 
(lower part of the Table) . Presumably, this affirmation is 
also valid for the electron-electron interaction of arbitrary 
range. The arguments justifying the downward-convexity 
condition for realistic quasi-one-dimensional compounds 
were presented by Hubbard ||l5| . 

Apart from the metal insulator transition, another in- 
teresting phenomenon studied in the spinless fermion 
model is the charge ordering. In particular, for p = 1/2 
the interplay of "usual" and dimerized (Wigner and 
Pcicrls |22|) charge ordering can be discussed. Our ap- 
proach cannot be applied directly to the solution of this 
problem since it is rather difficult to find asymptotics 
of the corresponding correlation functions directly from 
finite-cluster calculations. This problem can be inves- 
tigated either by such approaches as analytical or 
numerical ||25| renormalization group or by the combina- 
tion of exact diagonalization technique with the results 
of the theory of Luttinger liquid (see e.g. |^). Neverthe- 
less we present here our results about the characteristics 
of the short-range charge order in the model under con- 
sideration which can be interesting themselves. 

The results of the calculations for the correlation func- 
tions (noUi) and (710^2) in the ground state for p = 1/2 
are presented in Figs. 3 and 4. Consider first the case 
V < ^V. According to Fig. 3 in this case we always have 
(r^o?^l) < which show the absence of the dimeriza- 

tion. The metal-insulator transition is almost not appre- 
ciable in the calculated correlation functions describing 
the short-range order. 

Consider now the case V > ^V. The data presented 
in Fig. 4 show that as V is decreased, the dimerized 
state with {riQni) > (no"-2) is probably destroyed before 
the metal- insulator transition occurs. One may suppose 
that as the transfer integral t is increased starting from 
the atomic limit, first the dimer lattice melts and then 
the insulator-metal transition occurs. Unfortunately the 
present results give us only preliminary indications of the 
melting of dimerized lattice (because we cannot inves- 
tigate the asymptotics of the correlation functions and 
therefore have no direct information about long-range 
order), and this interesting question calls for further in- 
vestigations. 
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The results of calculations for the correlation functions 
(noni) , (710^2) and {non^} for p = 2/3 and = 10 is 
shown in Fig. 5. Qualitatively the similar picture takes 
place for any V < V. At large V and V (or equiv- 
alently for i — > 0) we have (noni) — > 1/3, (710/12) 
1/3, (nons) ^ 2/3. 

So our calculations give a rather full description of 
the metal-insulator transition in the ID spinless fermion 
model beyond the nearest neighbor approximation. They 
also provide a basis for further investigations in particular 
concerning the relation between the metal-insulator tran- 
sition and the destruction of the charge ordering. One 
may think that the derived conclusions about the effects 
of competing interactions on the nature of Verwey tran- 
sition may be useful in analyzing the realistic 2D or 3D 
systems. 



IV. ELECTRONIC TOPOLOGICAL TRANSITION 
IN THE MODEL WITH 
NEXT-NEAREST-NEIGHBOR TRANSFER 
PROCESSES 

Apart from the metal-insulator transition considered 
above the interelectron interaction may lead to the elec- 
tron phase transition of another type. It is connected 
with the next-nearest neighbor transfer processes. Al- 
though the role of such processes in many-electron mod- 
els is investigated in a number of works (see e.g. recent 
paper [2g|]) the opportunity of such transition, to our 
knowledge, was not noted and studied. 

Consider at first the case t' 0,V' = 0. In that case 
at V = one-particle electron spectrum is given by 



£fc ~ —2t cos k — 2t' cos 2k. 



(6) 



It is easily seen that at t' > t/A the spectrum is non- 
monotonic in the range (0,7r). At t' > t/2 and p — 1/2 
the second fermion "pocket" appears i.e. the electron- 
filled region in k space is no longer singly connected. In 
other words, while the Fermi surface of "normal" one- 
dimensional system consists of two points k = ikp it 
consists of four points at t' > t/2 and p ~ 1/2. In the 
latter case the occupied states corresponds to the inter- 
vals -TT < A: < ~k^, ~k^ <k<k% and fcj^ < fc < tt. It 



should be emphasized that ID gas of non-interacting elec- 
trons cannot exhibit a dispersion low of this type since 
it is well known (see e.g. |3^) that the one-dimensional 
one-electron Schroedinger equation for a periodic poten- 
tial must necessarily exhibit a monotone E{k) spectrum 
between k = and k = n (the lattice parameter is set 
equal to unity everywhere). Nevertheless, the Hamilto- 
nian (1) with t' > t/A may be used for the description of 
real systems e.g. consisting of the pairs of strongly cou- 
pled one-dimensional chains with weak coupling between 
different pairs. In this case t is the transfer integral be- 
tween the nearest sites in the direction across the double 
chain and t' is the transfer integral along it (see Fig. 6). 



Studying of such systems is not only of purely theo- 
retical interest. It may give a deeper insight into the 
properties of realistic quasi-one-dimensional systems in 
particular the well-known compound NMP-TCNQ (the 
molecules of TCNQ and NMP arc shown in Figs. 7a and 
7b, respectively). This compound is characterized by a 
charge transfer from TCNQ to NMP molecules (about 
1/3 electron per NMP molecule). As one NMP-TCNQ 
cell contains one electron, half of the k states between — tt 
and TT are occupied. Therefore the Fermi surface consists 
of two sheets: one of them is bounded by k^ — Tr/3 and 
corresponds to the TCNQ chain; the other is bounded by 
kp — 57r/6 and corresponds to NMP chain (see Fig. 7c) 

It is rather natural to discuss in such a case the fol- 
lowing question: What would happen to the distribution 
function upon switching the electron-electron interaction, 
i.e., would the increasing interaction result only in a grad- 
ual smearing of the steps in the distribution function of 
electrons in k space or a drastic change of its shape is also 
possible, e.g., the merging of two steps? It will be shown 
below that the second possibility really takes place. 

Before the presentation and discussion of exact numer- 
ical results it would be reasonable to treat this problem 
in the simplest Hartree-Fock approximation. The trial 
wave function will be chosen in the form corresponding 
to the distribution function with two steps in k space of 
the widths 2Trp — 2A and 2A: 



= 11410) 



(7) 



where quasimomentum takes on the values —irp+X < k < 
Trp—X (the first step) and —tt < k < — tt-I-A, tt— A < k < tt 
(the second step). 

Substituting (7) into (2) and replacing the sum- 
mation by the integration according to the formula 
(1/^) Efc ••■ ^ (l/27r) dk... one has 



E _ {i}\H\^) 
N ~ N 



p'^V t[sin{pTT - A) - sin A] 



- -t'[sm{2pTT - 2A) + sin2A] - -^[sin(p7r - A) - sin A]^. 

In the case of p = 1/2, the substitution x = 1 — cos A + 
sin A yields 



E V 2 2, ^, 

- = - - -t + -(i - 2i' + -)x + -(2t' 

N \ TT TT TT TT 



V 



)x\ (9) 



where 0<x<2,x = corresponds to one pocket and 
a; > corresponds to two pockets. Minimizing the en- 
ergy with respect to x one obtains that the increase of V 
results in the merging of the two steps (disappearing of 
the second pocket) at a critical value 



Vc = 1T{2t' -t). 



(10) 



This means that the switching on the electron-electron 
interaction, according to the usual Hartree-Fock approxi- 
mation, may lead to a topological transition i.e. a change 



(8) 
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in topology of Fermi surface: the two-sheeted (four point) 
"surface" transform into a one-sheeted (two-point) one. 
In a 3D case it would correspond to a transition from 
a doubly connected Fermi surface to a singly connected 
one which is a particular case of the electronic topological 
transitions proposed by Lifshitz jsj] . The investigation of 
such opportunity by a more rigorous way is the purpose 
of our calculations described below. 

The calculations have been carried out for a system of 
10 electrons in the ring from 20 sites and of 8 electrons 
in the ring from 16 sites. The results for the points of 
the electronic topological transition in these two sets of 
calculations were the same with the accuracy of two sig- 
nificant digits. The computational results are shown in 
Figs. 8,9. Fig. 8 displays the distribution function in the 
ground state (rifc) at t' = t. It is clearly seen that "the 
hump" near the edges of the Brillouin zone {k — — tt, tt) 
disappears at 2.5t < V < 2.6t. Fig. 9 presents the phase 
diagram of the system. The comparison with Eq.(lO) 
show that beginning from approximately t' w 0.8t the 
Hartree-Fock approximation gives too high values for Vc 
and the difference with numerical results grows with in- 
creasing t' . So, the Hartree-Fock approximation is not 
very accurate quantitatively at large V. 

In the case V ^ the inverse topological transi- 
tion turns out to be possible when the increase of the 
Coulomb interaction results in the appearance of the sec- 
ond pocket. We have carried out the corresponding cal- 
culations for the case t' = 0.49t, .p = 1/2. Thus, without 
interaction the second minimum in the electron spectrum 
lies very close to the Fermi level and above it. The in- 
teraction can make this minimum to be lower than the 
Fermi level. The phase diagram for this case is shown 
in Fig. 10. The Hartree-Fock approximation can describe 
qualiatively the phenomenon of " inverse" transition but 
gives too high values of the critical values of V (the dif- 
ference may be in a factor of order of 2-5 depending on 
V values). 

To conclude this section note why these results are, 
to our opinion, non-trivial. It is well-known that for ID 
systems the electron velocity may be zero either at /c = 
or at the boundary of the Brillouin zone |3^. There- 
fore, true one-dimensional systems cannot exhibit Van 
Hove singularities inside the allowed band and hence they 
exhibit no electronic topological transitions because the 
latter are nothing but the crossing of the Fermi level 
by the Van Hove singularity. However, this is not the 
case for, e.g., a double chain. The exact numerical 
results presented here suggest that in this case quasi- 
one-dimcnsional systems do may demonstrate electronic 
topological transitions. As in the 3D case |^ they are, 
unlike metal-insulator transitions, transitions of metal- 
metal type. Studying the possibility of such transitions 
in real quasi-one-dimensional conductors containing dou- 
ble chains would be of experimental interest. 

We are grateful to A.O.Anokhin and S.V.Tretjakov for 
the assistance in performing of numerical calculations. 
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CAPTIONS TO FIGURES 

Fig.l. Metal-insulator phase diagram at p = 1/2. 
Black and empty circles correspond to the metal and in- 
sulating phase, respectively; solid lines correspond to the 
condition Aqo = 2. 
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Fig. 2. Metal-insulator phase diagram at p = 2/3. 
Black and empty circles correspond to the metal and in- 
sulating phase, respectively; solid lines correspond to the 
condition Aqo = 3. 

Fig. 3. Dependences of the correlation functions (noni) 
(circles) and {11^712) (squares) on the parameter V for 
different V < V/2: (1) V = (2) V = 0.2V (3) V = 
0.3V (4) V = 0.5V. 

Fig. 4. Dependence of the correlation functions (noni) 
(circles) and {nQn2) (squares) on the parameter V for 

V = 0.7V. 

Fig. 5. Dependences of the correlation functions {riQni) 
(circles), (nori2) (squares) and {riQn^) (triangles) on V' 
at = 10t,p= 2/3. 

Fig. 6. The double chain. 

Fig. 7. (a) TCNQ molecule (hydrogen atoms are not 
shown) (b) NMP molecule (c) the sketch of the dispersion 
law of electrons in NMP-TCNQ according to . 

Fig. 8. The distribution function {ni^) at t' — t. The 
solid line corresponds to V — 0, the dash-dot line cor- 
responds to V = 2.5t, the dashed line corresponds to 

V = 2.6t. 

Fig. 9. Phase diagram of the ID spinless fermion model 
at p = 1/2, V' = 0, a being the label for "one-pocket" 
region and /3 for "two-pocket" one. 

Fig. 10. Phase diagram of the ID spinless fermion 
model a.t p — 1/2, t' — 0.49t, a being the label for "one- 
pocket" region and /3 for "two-pocket" one. 



TABLE I. Dependence of the ground state energy (the upper lines) and the gap A (the lower 
lines) on the different V and V at p = 2/3 



v 





1 


2 


3 


4 


5 


6 


8 


10 


V 



-0.5472 


-0.1318 


0.2662 


0.6489 


1.0188 


1.3790 


1.7324 


2.4248 


3.1072 




0.000 


-0.001 


-0.005 


0.004 


-0.004 


-0.003 


0.002 


0.004 


-0.001 


1 


-0.1856 


0.2347 


0.6386 


1.0276 


1.4031 


1.7675 


2.1236 


2.8180 


3.5016 




0.000 


0.000 


-0.003 


0.004 


0.003 


0.002 


0.001 


0.001 


-0.001 


2 


0.1669 


0.5903 


0.9984 


1.3922 


1.7728 


2.1417 


2.5010 


3.1997 


3.8839 




-0.001 


0.000 


0.000 


0.002 


-0.004 


0.002 


0.003 


0.003 


0.004 


3 


0.5135 


0.9395 


1.3495 


1.7459 


2.1292 


2.5008 


2.8626 


3.5662 


4.2545 




-0.002 


0.000 


0.000 


-0.001 


-0.004 


0.005 


0.005 


0.005 


0.006 


4 


0.8565 


1.2832 


1.6941 


2.0916 


2.4757 


2.8467 


3.2072 


3.9079 


4.5944 




-0.002 


0.000 


-0.003 


0.005 


0.094 


0.140 


0.169 


0.092 


0.083 


5 


1.1971 


1.6236 


2.0355 


2.4332 


2.8166 


3.1864 


3.5454 


4.2440 


4.9294 




-0.002 


0.000 


-0.001 


0.020 


0.309 


0.610 


0.876 


1.018 


0.932 


6 


1.5367 


1.9631 


2.3746 


2.7721 


3.1547 


3.5234 


3.8814 


4.5788 


5.2637 




-0.003 


0.000 


-0.003 


0.019 


0.374 


1.017 


1.497 


1.952 


1.934 


7 


1.8744 


2.3005 


2.7123 


3.1095 


3.4913 


3.8591 


4.2164 


4.9131 


5.5976 




-0.003 


0.000 


-0.004 


0.039 


0.422 


1.316 


2.042 


2.803 


2.925 


8 


2.2113 


2.6376 


3.0488 


3.4458 


3.8271 


4.1941 


4.5510 


5.2471 


5.9314 




-0.004 


0.000 


-0.001 


0.066 


0.586 


1.470 


2.447 


3.443 


3.822 


10 


2.8830 


3.3092 


3.7205 


4.1166 


4.4969 


4.8632 


5.2192 


5.9147 


6.5986 




-0.004 


0.000 


0.000 


0.082 


0.570 


1.390 


2.943 


4.568 


5.423 


20 


6.2307 


6.6539 


7.0639 


7.4582 


7.8362 


8.2005 


8.5555 


9.2496 


9.9330 




-0.004 


-0.003 


0.000 


0.127 


0.748 


1.793 


3.036 


5.781 


8.721 


100 


32.9074 


33.3307 


33.7389 


34.1312 


34.5074 


34.8702 


35.2243 


35.9175 


36.6004 




-0.005 


-0.003 


0.015 


0.173 


0.861 


1.963 


3.292 


5.983 


8.813 



TABLE IL Dependence of the gap A on the parameter A 
at different V, V and V" at p = 1/2 



V 


v 


V" 


Aoo 


A 


0.5 





0.25 


1.0 


0.000 


2.0 


1 


0.50 


1.0 


-0.001 


3.0 


2 


1.00 


1.0 


-0.001 


4.5 


3 


1.50 


1.5 


0.038 


1.0 





0.50 


2.0 


0.003 


4.0 


2 


1.00 


2.0 


0.096 


6.0 


4 


2.00 


2.0 


0.180 


7.5 


5 


2.50 


2.5 


0.465 


1.5 





0.75 


3.0 


0.263 


6.0 


3 


1.50 


3.0 


0.726 


9.0 


6 


3.00 


3.0 


0.873 


2.0 





1.00 


4.0 


1.069 


8.0 


4 


2.00 


4.0 


1.827 


10.0 


5 


2.50 


5.0 


3.119 


3.0 





1.50 


6.0 


3.335 


4.0 





2.00 


8.0 


5.296 


2.0 


2 


1.00 


2.0 


0.002 


3.0 


3 


1.50 


3.0 


0.000 


4.0 


4 


2.00 


4.0 


0.156 


5.0 


5 


2.50 


5.0 


1.769 


6.0 


6 


3.00 


6.0 


3.494 
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